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One dimensional chains of non-Abelian quasiparticles described by SU(2)k Chern-Simons-Witten 
theory can enter random singlet phases analogous to that of a random chain of ordinary spin- 1/2 
particles (corresponding to k — > oo). For k — 2 this phase provides a random singlet description of 
the infinite randomness fixed point of the critical transverse field Ising model. The entanglement 
entropy of a region of size L in these phases scales as Sl — log 2 L for large L, where d is the 
quantum dimension of the particles. 



A particularly exotic form of quantum order is possible 
in two space dimensions — so-called non-Abelian order 
In states with non-Abelian order, when certain local- 
ized quasiparticle excitations are present there is a low- 
energy Hilbert space whose dimensionality grows expo- 
nentially with the number of these quasiparticles. When 
these quasiparticles are well separated, this low-energy 
space becomes degenerate, and its states are character- 
ized by purely topological quantum numbers, meaning 
they cannot be distinguished by local measurements. If 
these quasiparticles are then adiabatically moved around 
one another, unitary transformations corresponding to 
non-Abelian representations of the braid group are car- 
ried out on this degenerate space. Aside from their in- 
trinsic scientific interest, recent attention has focused on 
the possibility of one day using non-Abelian states to 
perform fault-tolerant quantum computation 0, Q • 

Recently Feiguin et al. Q have studied models of 
interacting non-Abelian quasiparticles, specifically uni- 
form chains in which neighboring quasiparticles are close 
enough together to lift the degeneracy of the topologi- 
cal Hilbert space. In this Letter we study a related class 
of random interacting chains of non-Abelian quasiparti- 
cles. We are motivated both by [3| and by recent work of 
Refael and Moore [a, IfJ showing that the entanglement 
entropy of certain random one-dimensional models scales 
logarithmically with a universal coefficient. We find the 
same is true here for an infinite class of models. 

Exact diagonalization studies 0, @, 0] provide com- 
pelling evidence that the experimentally observed v = 
5/2 fractional quantum Hall (FQH) state is a non- 
Abelian state described by the Moore-Read "Pfaffian" 
state (H- This state belongs to a wider class of non- 
Abelian FQH states introduced by Read and Rezayi [l(| , 
labeled by index k. In this class, the k = 1 state is an 
ordinary (Abelian) Laughlin state, the k — 2 state is the 
Moore- Read state, and all subsequent integer fc values de- 
scribe new non-Abelian states. There is some numerical 
evidence 1(J U| that the k = 3 Read- Rezayi state de- 
scribes the experimentally observed v = 12/5 FQH state 
0. 

The quasiparticle excitations of the Read- Rezayi states 
with index k can be viewed (up to Abelian phases) as par- 
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FIG. 1: (a) Two non-crossing singlet states for SU(2)k par- 
ticles and their overlap, (b) Action of the singlet projection 
operators n? (which acts on particles 1 and 2) and II" (which 
acts on particles 2 and 3) on a particular non-crossing singlet 
state. The quantity d appearing in (a) and (b) is the quantum 
dimension of the particles. 



ticle excitations in SU(2)k Chern-Simons-Witten theory 
[l3j ]. These particles are characterized by their topolog- 
ical charge, a quantum number which can be viewed as 
a "<7-deformed" spin [14j. At level fc, topological charge 
can take the values 0,^,1, ■••,§, and obeys the fusion 
rule, 



Si <K> s 2 



min(si + s 2 , k - Si - s 2 ). (1) 



For k > 2 this implies \ <£> \ = © 1. Thus, when 
combining two particles with topological charge 1/2 the 
resulting state can either have topological charge or 
1. For ordinary spin-1/2 particles the former would be 
referred to as a singlet and the latter as a triplet. We will 
use the same terminology for SU{2)k particles, though 
it should be noted that here there is no S z degeneracy, 
i.e. there is only one "triplet" state. (For reviews of the 
general theory of non-Abelian particles see (Til fl6l|). 

The total spin sector of a one-dimensional chain of 
ordinary spin-1/2 particles is spanned by the the set of 
all "non-crossing" singlet states, i.e. states in which pairs 
of particles form singlet bonds in such a way that these 
bonds do not cross (see Fig. QJa)). Furthermore, these 



non-crossing states are linearly independent 17], and 
their number, and hence the dimensionality of the spin 
Hilbert space, grows asymptotically as 2 N for large N. 

Using the generalized notion of singlet described above, 
non-crossing singlet states can also be used as a basis for 
the total topological charge sector of a one-dimensional 
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FIG. 2: One step in the decimation procedure. 



chain of SU(2)k particles [18j. In this case the interpre- 
tation is that any pair of particles connected by a singlet 
bond will fuse to topological charge if brought together 

For N ordinary spin-1/2 particles, the overlap of two 
non-crossing singlet states \a) and \ft) can be computed 
by overlaying the two bond configuration and count- 
ing the number of closed loops, Ni oops . The overlap is 
then given by (a\/3) = 2 N '°°p°~ N / 2 . For SU{2) k parti- 
cles this overlap rule becomes (a\ft) — d Nlo °p s ~ N / 2 where 
d = 2 cos jrpj is the "quantum dimension" of the parti- 
cles (see Fig. Ufa)) [H For these values of d the 
non-crossing states are no longer linearly independent — 
they satis fy li near relations known as the Jones- Wenzl 
projectors [2(| which reduce the size of the Hilbert space 
so that its dimensionality grows asymptotically not as 
2 N , but as d N for large N. 

Consider a random one-dimensional chain of SU(2)k 
particles. Following Q, we assume that neighboring par- 
ticles are close enough together so that the singlet and 
triplet fusion channels are split in energy, with the sin- 
glet lying lowest. The Hamiltonian describing this chain 
is then 



(2) 



where Ji > is the energy splitting associated with par- 
ticles at sites i and i + 1, and 11° is the singlet projection 
operator on these particles, the action of which on repre- 
sentative non-crossing singlet states is shown in Fig.QJb). 
The uniform versions of these models (Ji — J) were stud- 
ied numerically for k — 3 and analytically for all k in [3| , 
where they were shown to be conformally invariant with 
central charge c = 1 — 6/((fc + l)(k + 2)). 

Because the Hilbert space of this SU(2)k chain can 
be described using a non-crossing singlet basis, the usual 
real-space renormalization group (RG) approach based 
on decimating singlet bonds [2l|, [22J can be straightfor- 
wardly applied to ([2|) when the J,*'s are random. Each 
iteration of this procedure begins by finding the strongest 
bond in the chain, i.e. the Ji with the highest value, and 
making the approximation that the two particles con- 
nected by it fuse to topological charge and so form a 
singlet bond. 

The effective interaction J between the two particles 
on either side of this singlet is then determined pertur- 
batively as follows. Consider four neighboring particles 
and the associated three bond strengths Ji, Ji and J3, 
with J 2 3> Ji, J3 so that, as described above, a singlet 



forms between the two particles connected by J2 (see 
Fig. [2]). A straightforward generalization of the usual 
second-order perturbation theory calculation for ordinary 
spin-1/2 particles, but using the modified overlap rules 
shown in Fig. [TJ then yields, 



J = (2/d 2 )J 1 J 3 /J 2 . 



(3) 



Provided d > y/2, which is the case for all k > 2 consid- 
ered here, J will always be less than the strength of the 
decimated bond J 2 . Thus, as this procedure is iterated, 
high-energy bonds are systematically eliminated, leading 
eventually to a single non-crossing singlet state. 

The RG flow produced by this decimatio n p rocedure 
can then be analyzed in the standard way 2l|, |22| . In- 
troducing the logarithmic bond strength variables Pi = 
ln(fi/Ji), where ft is the largest remaining bond strength 
at any given stage of decimation, ([3]) can be written 
ft = ft\ + /?3 + \n(2/d 2 ). Defining the flow parameter 
r = ln(f2o/f2) where fio is the largest bond strength 
at the outset of the decimation procedure, and ignor- 
ing the ln(2/d 2 ) (i.e. taking ft ~ ftx + fta), an approx- 
imation which can be justified a posteriori due to the 
broad distribution of /3's at the fixed-point, an integro- 
differential equation can be written down for the bond 
strength distribution Pr (ft) [U • This distribution is 
defined so that when the flow parameter is T the frac- 
tion of bonds with logarithmic strength between ft and 
ft + dft is Pv(ft)dft. As shown by Fisher 0, for almost 
any initial random bond configuration, the bond strength 
distribution flows to the infinite randomness fixed point 
distribution, Pv(ft) = e~^/ r /Y. The resulting phase is 
known as a random singlet phase. 

It follows that the random SU(2)k chains |2]) flow to 
random singlet phases for all k > 2. In the limit fc — > 00 
this phase corresponds to the usual random singlet phase 
for ordinary spin-1/2 particles [22T |. For k = 2 we now 
show that the resulting phase can be mapped onto the 
infinite randomness fixed point of the critical transverse 
field Ising model [23)], thus providing a "random singlet" 
description of this fixed point. 

We use the fact that 5 , C/(2) 2 particles can be repre- 
sented using Major ana fermions operators 74 [l(| — op- 
erators which are self conjugate (jj — 7$) and which 
satisfy the Clifford algebra {ji,jj} — 2§ij. Two 
Majorana fermions can be combined to form a usual 
fermion, so that, e.g., associated with neighboring sites 



and i + 1 there is a fermion operator c. 



in 



z7i+i)/-\/2 which satisfies the usual anticommutation re- 
lation {ci t i+i,c\ — 1 and which anticommutes with 
any similar fermion operator constructed out of a differ- 
ent pair of Majorana fermions. The Fermi mode associ- 
ated with this pair can then be occupied (corresponding 
to topological charge 1) or unoccupied (corresponding to 
topological charge 0). The singlet projection operator is 
then IT? = 1 — c\ i+1 Ci^ + i, which in the Majorana repre- 
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FIG. 3: "Random singlet" view of a decimated random trans- 
verse field Ising model. A random singlet state (green) is 
overlaid with a dimer state (blue). In the dimer state bonds 
connect pairs of SU(2)2 particles which are mapped onto the 
spin- 1/2 degrees of freedom of the transverse field Ising model 
(bottom row of dots). Closed loops then correspond to deci- 
mated superspins, indicated by solid arrows. 




FIG. 4: Schmidt decomposition of a state of N pairs of SU (2)h 
particles connected by singlet bonds. The states in the decom- 
position are expressed using a basis in which circles enclose 
particles in topological charge eigenstates. The sum is over 
all S2, S3 • • • , sjv consistent with the fusion rule (fl]). 



sentation is IT? = 

To map the SU (2)2 chain onto the transverse field Ising 
model we first group together neighboring pairs of Majo- 
rana fermions. Letting the index j label these pairs, each 
of which consists of a right Majorana fermion (7?) and a 
left Majorana fermion (7^), the Hamiltonian ([2]) can be 
written 



H 



E» 



E 
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(4) 



Here hj corresponds to the coupling within the j pair, 
and Jj corresponds to the coupling between the right- 
most particle in the j th pair and the leftmost particle 
in the (j + l) st pair. The usual Jordan- Wigner trans- 
formation (see, for example, ) , jj = erj IIfc=i a k an< 4 

Tf = °j ni=i °fc ! tnen ma PS © onto the random trans- 
verse field Ising model, 



H E h rf • E 



(5) 



Because hj and Jj are drawn from the same distribution, 
the model is at its critical point. 

The usual decimation procedure for the transverse field 
Ising model involves two separate steps — either forming 
ever larger "superspins" when the strongest interaction 
is an Ising interaction (J), or decimating these super- 
spins when the strongest interaction is a magnetic field 
strength (h) (H- Tnc SU(2) 2 "random singlet" view of 
this decimation provides a unified description of these 
two steps. Figure [3] shows a random singlet state pro- 
duced by decimation and a reference "dimer" state in 
which bonds connect pairs of Majorana fermions which 
correspond to single spins in the transverse field Ising 
model. Overlaying these two states produces closed loops 
which, in the transverse field Ising model, correspond 
to decimated superspins. Essentially, as the decimation 
which produces the random singlet state shown in the 
figure proceeds, any time a bond forms which does not 
close a loop this corresponds to eliminating an Ising in- 
teraction and increasing the number of spins contributing 



to a superspin. Then, when a bond forms which closes 
a loop, the corresponding superspin is frozen along the 
direction of the applied field and decimated. 

Recently Refael and Moore [B| have shown that the en- 
tanglement entropy associated with the infinite random- 
ness fixed points of both the random spin-1/2 Heisen- 
berg chain (fc — > 00) and the transverse field Ising model 
(fc = 2) have universal scaling properties which can be 
used to generalize the notion of central charge to one- 
dimensional quantum critical systems which are not con- 
formally invariant. We now show that the same is true for 
all the SU(2)k infinite randomness fixed points. The en- 
tanglement entropy of these states is calculated by treat- 
ing a contiguous segment of the chain consisting of L 
particles as a subsystem (denoted A) of the full chain. 
Tracing out the degrees of freedom of the rest of the chain 
then yields a reduced density matrix pa- The entangle- 
ment entropy is the average over realizations of disorder 
of the Von Neumann entropy of this reduced density ma- 
trix, S L = - (Tr p A log 2 pa) ■ 

In random singlet states the calculation of Sl for large 
L can be done, as in by counting the number of singlet 
bonds which connect sites in region A with sites outside 
of it, averaging over realizations of disorder, and then 
multiplying the result by the entanglement entropy asso- 
ciated with each bond. All the SU(2)k random singlet 
states discussed here are governed by the same fixed point 
bond distribution as that considered in so the result 
of that work that the average number of bonds contribut- 
ing to the entanglement scales as | In L for large L holds 
here as well. 

To compute the entanglement entropy per bond for 
SU(2)k particles, imagine forming N singlet pairs, with 
one particle from each pair taken to be in subsystem A, 
the other in subsystem B, as shown in Fig. [5] This fig- 
ure also shows a Schmidt decomposition of this state us- 
ing a basis in which ovals are drawn around particles in 
topological charge eigenstates. The Schmidt coefficients 
(A Sjv in Fig. [4| can be obtained using standard calcula- 
tion techniques for non-Abelian particles 15, lj|. They 
depend only on the total topological charge sjv of the 
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particles in region A (or equivalently region B) of the 
corresponding state in the Schmidt decomposition, and 
are given by X SN — [2s n + l]/d N , where we have intro- 
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of the reduced den- 
tracing out the de- 
is then El, S A = 



duccd the g-integers 
with q = cxp il-K l(k + 2). 

The Von Neumann entropy 
sity matrix pa obtained by 
grees of freedom in region B 
- J2 SN s N )X SN log 2 X SN , where D(JV, sat) is the di- 
mensionality of the space of N SU(2)k particles with to- 
tal topological charge sn- Using the fact that, for large 
NJXN,s n ) ~ [2s N +l]d N /V 2 where V 2 = £^ 2 [2s+l] 2 
[IE El, it follows that S A - Nlog 2 d- 0(log 2 k) for 
N 3> k. Thus for large N the entanglement per bond 
is log 2 d, reflecting the fact that the size of the Hilbert 
space of N particles grows asymptotically as d N [261 ]. 

Returning to the SU{2)k random singlet phases, mul- 
tiplying the average number of bonds leaving a region of 
size L (~ i In L) by the entanglement per bond (~ log 2 d) 
yields 



S L ~(lnd/3)log 2 L. 



(6) 



Following [g], if we compare (O with the entangle- 
ment entropy of conformally invariant one-dimensional 
systems. S 



| log 2 L where c is the central charge 
it is natural to define an "effective cen- 
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tral charge" of c = \nd for these phases. In the k — ► oo 
limit, corresponding to the ordinary SU{2) random sin- 
glet phase with d = 2, we have c = In 2, and for k = 2, 
corresponding, as shown above, to the critical point of 
the random transverse field Ising model with d = y2, we 
have c = 5 In 2, both of which agree with results obtained 
in@. 

Finally we note that for the SU(2)k chains considered 
here the effective central charge of the disordered model, 
c = Ind, is always less than the central charge of the 
uniform model, c = 1 - 6/((fc + l)(fc + 2)) [J, though 
the simple relation c = In 2 x c emphasized in Q only 
holds for k — > 00 and k = 2. This is consistent with the 
generalized "c-theorem" envisioned in [B[ which supposes 
that the effective central charge decreases along RG flows 
between quantum critical points. However, it should be 
emphasized that this "theorem" is not a rigorous result. 
In particular, Santachiara |30| has shown that it is vio- 
lated by RG flows from the uniform to disordered phases 
of the Z n parafermionic Potts model for n > 42. 
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